Explicit modeling of volatile organic compounds partitioning in the atmospheric aqueous phase

The gas phase oxidation of organic species is a multigenerational process involving a large number of secondary compounds. Most secondary organic species are water-soluble multifunctional oxygenated molecules. The fully explicit chemical mechanism GECKO-A (Generator of Explicit Chemistry and Kinetics of Organics in the Atmosphere) is used to describe the oxidation of organics in the gas phase and their mass transfer to the aqueous phase. The oxidation of three hydrocarbons of atmospheric interest (isoprene, octane and α-pinene) is investigated for various NOx conditions. The simulated oxidative trajectories are examined in a new two dimensional space defined by the mean oxidation state and the solubility. The amount of dissolved organic matter was found to be very low (yield less than 2% on carbon atom basis) under a water content typical of deliquescent aerosols. For cloud water content, 50% (isoprene oxidation) to 70% (octane oxidation) of the carbon atoms are found in the aqueous phase after the removal of the parent hydrocarbons for low NOx conditions. For high NOx conditions, this ratio is only 5% in the isoprene oxidation case, but remains large for α-pinene and octane oxidation cases (40% and 60%, respectively). Although the model does not yet include chemical reactions in the aqueous phase, much of this dissolved organic matter should be processed in cloud drops and modify both oxidation rates and the speciation of organic species

Modeling SOA formation from the oxidation of intermediate volatility <i>n</i>-alkanes

The chemical mechanism leading to SOA formation and ageing is expected to be a multigenerational process, i.e. a successive formation of organic compounds with higher oxidation degree and lower vapor pressure. This process is here investigated with the explicit oxidation model GECKO-A (Generator of Explicit Chemistry and Kinetics of Organics in the Atmosphere). Gas phase oxidation schemes are generated for the C₈–C₂₄ series of <i>n</i>-alkanes. Simulations are conducted to explore the time evolution of organic compounds and the behavior of secondary organic aerosol (SOA) formation for various preexisting organic aerosol concentration (<i>C</i>_OA). As expected, simulation results show that (i) SOA yield increases with the carbon chain length of the parent hydrocarbon, (ii) SOA yield decreases with decreasing <i>C</i>_OA, (iii) SOA production rates increase with increasing <i>C</i>_OA and (iv) the number of oxidation steps (i.e. generations) needed to describe SOA formation and evolution grows when <i>C</i>_OA decreases. The simulated oxidative trajectories are examined in a two dimensional space defined by the mean carbon oxidation state and the volatility. Most SOA contributors are not oxidized enough to be categorized as highly oxygenated organic aerosols (OOA) but reduced enough to be categorized as hydrocarbon like organic aerosols (HOA), suggesting that OOA may underestimate SOA. Results show that the model is unable to produce highly oxygenated aerosols (OOA) with large yields. The limitations of the model are discussed

Explicit modelling of SOA formation from α-pinene photooxidation: sensitivity to vapour pressure estimation

The sensitivity of the formation of secondary organic aerosol (SOA) to the estimated vapour pressures of the condensable oxidation products is explored. A highly detailed reaction scheme was generated for α-pinene photooxidation using the Generator for Explicit Chemistry and Kinetics of Organics in the Atmosphere (GECKO-A). Vapour pressures (P^(vap)) were estimated with three commonly used structure activity relationships. The values of P^(vap) were compared for the set of secondary species generated by GECKO-A to describe α-pinene oxidation. Discrepancies in the predicted vapour pressures were found to increase with the number of functional groups borne by the species. For semi-volatile organic compounds (i.e. organic species of interest for SOA formation), differences in the predicted Pvap range between a factor of 5 to 200 on average. The simulated SOA concentrations were compared to SOA observations in the Caltech chamber during three experiments performed under a range of NO_x conditions. While the model captures the qualitative features of SOA formation for the chamber experiments, SOA concentrations are systematically overestimated. For the conditions simulated, the modelled SOA speciation appears to be rather insensitive to the P^vap estimation method

Explicit modeling of volatile organic compounds partitioning in the atmospheric aqueous phase

The gas phase oxidation of organic species is a multigenerational process involving a large number of secondary compounds. Most secondary organic species are water-soluble multifunctional oxygenated molecules. The fully explicit chemical mechanism GECKO-A (Generator of Explicit Chemistry and Kinetics of Organics in the Atmosphere) is used to describe the oxidation of organics in the gas phase and their mass transfer to the aqueous phase. The oxidation of three hydrocarbons of atmospheric interest (isoprene, octane and α-pinene) is investigated for various NOx conditions. The simulated oxidative trajectories are examined in a new two dimensional space defined by the mean oxidation state and the solubility. The amount of dissolved organic matter was found to be very low (yield less than 2% on carbon atom basis) under a water content typical of deliquescent aerosols. For cloud water content, 50% (isoprene oxidation) to 70% (octane oxidation) of the carbon atoms are found in the aqueous phase after the removal of the parent hydrocarbons for low NOx conditions. For high NOx conditions, this ratio is only 5% in the isoprene oxidation case, but remains large for α-pinene and octane oxidation cases (40% and 60%, respectively). Although the model does not yet include chemical reactions in the aqueous phase, much of this dissolved organic matter should be processed in cloud drops and modify both oxidation rates and the speciation of organic species

Estimation of rate coefficients for the reactions of O3 with unsaturated organic compounds for use in automated mechanism construction

Reaction with ozone (O3) is an important removal process for unsaturated volatile organic compounds (VOCs) in the atmosphere. Rate coefficients for reactions of 3with VOCs are therefore essential parameters for chemical mechanisms used in chemistry transport models. Updated and extended structure-activity relationship (SAR) methods are presented for the reactions of O3 with mono- and polyunsaturated organic compounds. The methods are optimized using a preferred set of data including reactions of O3 with 221 unsaturated compounds. For conjugated dialkene structures, site-specific rates are defined, and for isolated polyalkenes rates are defined for each double bond to determine the branching ratios for primary ozonide formation. The information can therefore guide the representation of the O3 reactions in the next generation of explicit detailed chemical mechanisms

Estimation of rate coefficients and branching ratios for gas-phase reactions of OH with aromatic organic compounds for use in automated mechanism construction

Reaction with the hydroxyl (OH) radical is the dominant removal process for volatile organic compounds (VOCs) in the atmosphere. Rate coefficients for the reactions of OH with VOCs are therefore essential parameters for chemical mechanisms used in chemistry transport models, and are required more generally for impact assessments involving estimation of atmospheric lifetimes or oxidation rates for VOCs. A structure–activity relationship (SAR) method is presented for the reactions of OH with aromatic organic compounds, with the reactions of aliphatic organic compounds considered in the preceding companion paper. The SAR is optimized using a preferred set of data including reactions of OH with 67 monocyclic aromatic hydrocarbons and oxygenated organic compounds. In each case, the rate coefficient is defined in terms of a summation of partial rate coefficients for H abstraction or OH addition at each relevant site in the given organic compound, so that the attack distribution is defined. The SAR can therefore guide the representation of the OH reactions in the next generation of explicit detailed chemical mechanisms. Rules governing the representation of the reactions of the product radicals under tropospheric conditions are also summarized, specifically the rapid reaction sequences initiated by their reactions with O2

Estimation of mechanistic parameters in the gas-phase reactions of ozone with alkenes for use in automated mechanism construction

Reaction with ozone is an important atmospheric removal process for alkenes. The ozonolysis reaction produces carbonyls and carbonyl oxides (Criegee intermediates, CI), which can rapidly decompose to yield a range of closed shell and radical products, including OH radicals. Consequently, it is essential to accurately represent the complex chemistry of Criegee intermediates in atmospheric models in order to fully understand the impact of alkene ozonolysis on atmospheric composition. A mechanism construction protocol is presented which is suitable for use in automatic mechanism generation. The protocol defines the critical parameters for describing the chemistry following the initial reaction, namely the primary carbonyl/CI yields from the primary ozonide fragmentation, the amount of stabilisation of the excited CI, the unimolecular decomposition pathways, rates and products of the CI, and the bimolecular rates and products of atmospherically important reactions of the stabilised CI (SCI). This analysis implicitly predicts the yield of OH from the alkene–ozone reaction. A comprehensive database of experimental OH, SCI and carbonyl yields has been collated using reported values in the literature and used to assess the reliability of the protocol. The protocol provides estimates of OH, SCI and carbonyl yields with root mean square errors of 0.13 and 0.12 and 0.14, respectively. Areas where new experimental and theoretical data would improve the protocol and its assessment are identified and discussed

Exploration of the influence of environmental conditions on secondary organic aerosol formation and organic species properties using explicit simulations: development of the VBS-GECKO parameterization

Atmospheric chambers have been widely used to study secondary organic aerosol (SOA) properties and formation from various precursors under different controlled environmental conditions and to develop parameterization to represent SOA formation in chemical transport models (CTMs). Chamber experiments are however limited in number, performed under conditions that differ from the atmosphere and can be subject to potential artefacts from chamber walls. Here, the Generator for Explicit Chemistry and Kinetics of Organics in the Atmosphere (GECKO-A) modelling tool has been used in a box model under various environmental conditions to (i) explore the sensitivity of SOA formation and properties to changes on physical and chemical conditions and (ii) develop a volatility basis set (VBS)-type parameterization. The set of parent hydrocarbons includes n-alkanes and 1-alkenes with 10, 14, 18, 22 and 26 carbon atoms, α-pinene, β-pinene and limonene, benzene, toluene, o-xylene, m-xylene and p-xylene. Simulated SOA yields and their dependences on the precursor structure, organic aerosol load, temperature and NOx levels are consistent with the literature. GECKO-A was used to explore the distribution of molar mass, vaporization enthalpy, OH reaction rate and Henry's law coefficient of the millions of secondary organic compounds formed during the oxidation of the different precursors and under various conditions. From these explicit simulations, a VBS-GECKO parameterization designed to be implemented in 3-D air quality models has been tuned to represent SOA formation from the 18 precursors using GECKO-A as a reference. In evaluating the ability of VBS-GECKO to capture the temporal evolution of SOA mass, the mean relative error is less than 20&thinsp;% compared to GECKO-A. The optimization procedure has been automated to facilitate the update of the VBS-GECKO on the basis of the future GECKO-A versions, its extension to other precursors and/or its modification to carry additional information.</p

Structure-activity relationships to estimate the effective Henry's law constants of organics of atmospheric interest

International audienceAbstract. The Henry's law constant is a key property needed to address the multiphase behaviour of organics in the atmosphere. Methods that can reliably predict the values for the vast number of organic compounds of atmospheric interest are therefore required. The effective Henry's law constant H* in air-water systems at 298 K was compiled from literature for 488 organic compounds bearing functional groups of atmospheric relevance. This data set was used to assess the reliability of the HENRYWIN bond contribution method and the SPARC approach for the determination of H*. Moreover, this data set was used to develop GROMHE, a new Structure Activity Relationship (SAR) based on a group contribution approach. These methods estimate logH* with a Root Mean Square Error (RMSE) of 0.38, 0.61, and 0.73 log units for GROMHE, SPARC and HENRYWIN respectively. The results show that for all these methods the reliability of the estimates decreases with increasing solubility. The main differences among these methods lie in H* prediction for compounds with H* greater than 103 M atm−1. For these compounds, the predicted values of logH* using GROMHE are more accurate (RMSE = 0.53) than the estimates from SPARC or HENRYWIN